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ABSTRACT:  The  minimum  root  separation  of  a  complex  polynomial  A  is 
defined  as  the  minimum  of  the  distances  between  distinct 
roots  of  A.  For  polynomials  with  Gaussian  integer  coeffi¬ 
cients  and  no  multiple  roots,  three  lower  bounds  are  de¬ 
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absolute  values  of  the  coefficients  of  A.  The  notion  of 
a  semi-norm  for  a  commutative  ring  is  defined,  and  it  is 
shown  how  any  semi-norm  can  be  extended  to  polynomial 
rings  and  matrix  rings,  obtaining  a  very  general  analogue 
of  Hadamard's  determinant  theorem. 
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THE  MINIMUM  ROOT  SEPARATION  OF  A  POLYNOMIAL 


1.  Introduction.  Let  A(x)  be  a  polynomial  of  degree  nX  with 
complex  coefficients  a^  and  complex  roots  <y  >  so  that 


A(x)=^  "=0  aixl=anTTj  =  l(x-°j)- 


(0 


We  define  sep(A)  ,  the  minimum  root.  separation  of  A,  by 


sep(A'=min  ,  U.-aul«  ‘  ’ ' 

O'jfO'k  J  K 

with  the  convention  that  sep(A)=  ®  in  case  A  lias  only  one  distinct  root. 

The  computing  time  required  by  any  algorithm  to  isolate  the  zeros 
of  A  depends  inversely  on  sep(A)  .  Hence  we  are  interested  in  easily  com¬ 
putable  functions  f(a  ,...,an>  of  the  coefficients  such  that 

0<f(aQ, ...,an)<  sep'A)  .  ('0 

Heindel,  [ j] ,  in  analyzing  the  computing  time  of  an  algorithm 
based  on  Sturm's  theorem  for  isolating  the  real  zeros  of  any  polynomial 
with  integer  coefficients,  used  a  weak  lower  bound  for  sep(A)  due  to 
Collins.  Pinkert,  [9],  presents  an  analogous  algorithm  for  isolating  all 
zeros,  real  and  complex,  of  any  polynomial  with  Gaussian  integer  coefficients. 
His  algorithm  is  based  on  Sturm's  theorem  and  the  Routh-Hurwitz  theorem  and 
uses  a  stronger  lower  bound  for  sep(A)  obtained  more  recently  by  Collins. 
Horowitz,  using  another  simpler  approach,  has  recently  obtained  a  third 
lower  bound,  intermediate  in  strength,  but  just  slightly  weaker  than  the 


stronger  bound  of  Collins.  In  the  following,  these  three  bounds  are  all 
derived,  with  the  hope  of  stimulating  further  research  on  the  problem. 

If  A(x)  has  rational  comolex  coefficients,  we  can  easily  compute 
another  polynomial,  having  the  same  roots,  with  Gaussian  integer  co¬ 
efficients.  Further,  if  A(x)  has  Gaussian  integer  coefficients,  we  can 
easily  compute  another  polynomial  A*(x)  with  Gaussian  integer  coefficients, 
having  the  same  roots  as  A(x)  and  having  only  simple  roots,  namely 

A*(x)-A(x)/gcd(A(x),A'(x)),  (M 

where  A'(x)  is  the  derivative  of  A(x)  '‘nd  "ged"  denotes  the  greatest  common 
divisor.  Hence  in  the  following  A  is  assumed  to  haoe  GaLSSian  integer  co¬ 
efficients  arid  no  multiple  roots. 

Also,  the  three  lower  bounds  to  be  obtained  will  all  be  of  the  forr 

0<g(n,d)£sep(A),  (5) 

where  n=deg(A) ,  the  degree  of  A,  and  d=v(A) ,  where  v  is  some  "semi-norm". 

In  the  next  section  we  introduce  the  notion  of  a  semi-norm  for  a  ring  and 
then  derive  some  lemmas  which  will  be  used  in  deriving  the  root  separation 
theorems . 


2.  Semi-Norms  and  Resultants.  If  $  is  any  commutative  ring,  a 
semi-norm  for  p  is  any  function  v  from  into  the  non-negative  real  numbers 
satisfying  th<  following  three  conditions  for  all  a,  b^. 


\)( a) -0  if  and  only  if  a=0,  (f  a 

v(a-b^<  \/a )'*v(b ) »  0  b 

vf  ab  <  v(a )  v(  d)  •  0  c ) 

These  conditions  imply  also 

v(  -a  ,=\j/a)  ,  (6d 

v(a+b  _  \)\  a  +v(b).  O'0; 


A  norm  for  f?  is  a  semi-norm  for  ft  such  that 


v (ab  *vfa)v(b) 


C7) 


For  the  ring  G  of  the  Gaussian  integers  a  fanuliar  norm  is 

o  1  jo 

v(a+bi)s=(a+bil!=,a‘ +b"‘)  .  A  semi-norm  for  G  which  is  not  a  norm  is 

■  ( a+bi  )*=  |  a+bi  |  =  |  a  j  +  |  b  j  . 

Any  semi-norm  v  on  a  commutative  ring  can  be  extended  to  a  semi¬ 
norm  on  the  polynomial  ring  $x]  by  the  definition 

i=OaiX  ='  i=0v  ai'  ^  ' 

By  induction  on  r,  repeated  application  of  (£)  extends  v  to  a  semi-norm  on 
/^fx^,...,x  ],  which  is  easily  seen  to  be  independent  of  the  order  in  which 

the  indeterminates  x^  are  adjoined. 

As  a  special  case,  defines  j  A  I  and  |a|^  for  any  Gaussian  p  o ly ~ 

nomial  A(x1? . . .  ,xr)  eGfx^  . . .  .x^]  as  extensions  of  the  semi-norms  for  G 
defined  above.  For  integral  polynomials  Avx^,...,x^  with  rational  integer 
coefficients,  the  norm  |A|^  has  been  used  extensively  for  the  analysis  of 


I 


algebraic  algorithms.  See,  for  example,  [ 1] ,  [2  ,  [~  and  [6].  Its 
extension  to  Gaussian  polynomials,  however,  is  new. 

If  M  is  an  arbitrary  matrix  (or  vector  over  ft,  we  define 

v(M)^i2]v(Mij},  (9) 

where  the  summation  extends  over  all  entries  of  M.  It  is  easy  to  verify 
that  the  conditions  (6a)-(6c)  hold  tor  matrices  over  /^whenever  the  opera¬ 
tions  are  defined.  In  pa-ticular,  this  extends  v  to  a  semi-norm  for  the 
ring  of  all  n  by  n  square  matrices  over  ft. 

By  combining  the  semi-norm  extensions  for  polynomials  and  matrices, 
we  obtain  the  following  general  analogue  of  Hadamard 's  determinant  theorem 
([6],  p.  208). 


Theorem  1.  Let  ft  be  a  commutative  ring,  v  a  semi-norm  for  ft,  M  an 
n  by  n  matrix  over  ft.  Then 

v(det^M))<  TT"=1v(M.)  (10) 

where  M.  is  the  i^  row  of  M  and  det(M)  is  the  determinant  of  M. 

Proof .  By  induction  on  n,  the  case  n=l  being  trivial.  We  denote 
by  M.  .  the  element  of  M  in  the  i^  row  and  column  of  M,  by  M1.  .  the 

y  i,J  i.J 

submatrix  of  M  obtained  by  deletion  of  the  i*”*1  row  and  j1"  column.  By 
Laplace  expansion, 

det(M)  =  T  ^J(-l)J+1M1.det(M’1.).  (11) 

By  (6)  and  (11), 

v(det'M))<  ^  j^|v(M1j)v(det(M,1j))  . 


(12) 
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The  ith  rnw  of  M'^  is  a  subrow  of  Mi+1,  so 

v(det(M 1 1  j )  )<  TtJ^v(M1) 

U5) 

by  the  induction  hypothesis.  By  (12)  and  (13) » 

«(dot  (M))<  <>!&<%>• 

(!’•) 

Since  ^|v(Ml  .  )=v(Mj_) ,  this  completes  the  induction.! 

A  corollary  of  Theorem  1,  needed  in  Section  3,  will  now  be  obtained  by 
consideration  of  certain  submatrices  of  the  Sylvester  matrix  of  two  poly¬ 
nomials,  A  and  B,  over  Let  m=deg(A),  n=deg  B  .  The  Sylvester  matrix 
of  A  ari  B  is  the  m+n  by  m+n  matrix  S  whose  successive  rows  are  the  co¬ 
efficients  of  the  polynomials  xP  ^A(x) . xA(x) ,  A(x),X  B,x’,  ...» 

xB(x),  B(x)  .  Diagrammatically,  if  A(x)=^=(?aixi  and  B  x)=j=0bix  ,  then 


in  which  all  missing  entries  are  zero.  By  definition,  the  resultant  of 
A  and  B,  res(A.B),  is  the  determinant  of  S. 

Theorem  2.  Let  A  and  B  be  polynomials  over  a  commutative  ring  p 
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with  semi-norm  v  Let  m=degfA)>0,  n-deg;B)>0,  c-res  A,B  Then 

v(c)<v(A)  v(B)m. 

Also,  there  exist  polynomials  U  and  V  over  such  that  AUrBV=c,  deg(U)<n, 
deg(V)<m, 

v(u)<v(A)n_1v(B)m,  (17) 

and 

v(V)<v(A)nv(B)m_1.  (16) 

A.  I- 

Proof.  If  S.  is  the  i  row  of  S  then  v(S.)=v(A)  for  l<i<n  and 

-  L  1  ~ 

v(S^)=v(B)  for  n+l<i<  m+n,  and  (16)  follows  from  Theorem  1.  Now  consider 

y  m+n- i 

the  matrix  S  which  is  obtained  by  adding  to  the  last  column  of  S  x 

times  the  i^column  Qf  S,  for  l<i<m+n.  det(S  )=det(S)=c  and  the  last 

■X'  n — 1  /  \  <  \  *  /  \  ni“l  /  \ 

column  of  S  contains  the  entiles  x  A(x) , . . .  ,xA(x/>  A(x;,  x  B(xj . 

xB(x),B(x).  Applying  the  Laplace  determinant  expansion  theorem  to  the 

last  column  of  S*  we  obtain  AU+BV=e  with  deg(U)^n-l  and  deg(V)<m-l,  where 

x 

the  coefficients  of  U  and  V  are  the  cofactors  of  the  last  column  of  S  . 
Each  coefficient  of  U  is  the  determinant  of  a  matrix  obtained  from  E 
by  deleting  one  row  of  coefficients  of  A  and  the  last  column,  and  so 
Theorem  1  yields  (17),  and  similarly  (18)  holds.! 

5.  Root  Separation  Bounds.  For  each  of  the  first  two  root 
separation  bounds  we  will  i.se  the  following  upper  bound  on  the  roots  of 
a  polynomial. 

Theorem  Let  A  be  any  non-zero  Gaussian  polynomial,  and  let  c£, 
be  a  root  of  A.  Then 

W<|Al/laJ 

where  1  =ldcf'  A^  . 
n  ' 


(19) 


Proof .  |Al>la  I,  so  (19)  holds  for  |^|<1.  Let  A  x)~^_  a^x  and 

■  n  ©n-l  i 

assume  lry|^l.  Then  a  a  =l_j_=Qaia  »  so 

Dividing  (20)  by  \  ^ I n  ^ , 


(20) 


from  which  (19)  is  immediate.! 

Theorem  h .  (Collins,  I97O)  Let  A  be  a  Gaussian  polynomial  of 
degree  n^2  with  only  simple  roots,  and  d= | A I  .  Then 

sep(A)>(2d)"n(n'lW°.  (22) 


Proof.  Let  q,  ,.  ...^  be  the  zeros  of  A  and  A=sep(A).  We  may 
choose  notation  so  that  A^q^t^I*  Let  D  he  ^he  discriminant  of  A, 
so  that 


„  2ri-2 

D=a  T T. 


j<k 


r 

ai'ak^ 


(23) 


and  ([10],  Section  28),  D  is  a  Gaussian  integer.  Since  the  cxj  are 
distinct,  Da©  and  hence  |d|&1.  Combining  this  with  (23),  we  have 

<2"> 

Dividing  by  A  , 

Iv^'2-  (25) 

(J,k)V(l,2) 

There  are  (n  -n-2)/2  factors  loj_0'kl^  in  (25)  and  I  Oj  |  +  | 

2d/|a^|  by  Theorem  3*  Hence, 

2  2 
2  /r,j\n  -n-2  / 1  |H  ~3n 

A  *(2d)  /lan' 

2 

Now  n  -3n+2^0  and  la  |?1  so 


(26) 
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A'2s(2d)n2-n'2lanl2<(2d)n2*n, 


(27) 


from  which  (22)  is  immediate.  I 

Theorem  (Horowitz,  1973)  Let  A  be  a  Gaussian  polynomial  of 

degree  n^2  with  only  simple  roots,  and  d= | A |  .  Then 

sep(A)a(nd)  ^ .  (28) 

Proof.  Le t  .  ,oCn  oe  the  zeros  of  A  and^\=sep(A).  We  may 

suppose  that  X  ~ JL2 1 •  ®y  Theorem  2,  there  exist  Gaussian  poly¬ 
nomials  U  and  V  such  that 

AU+A'V=c,  (29) 

deg(u)<Jl-2  and  deg(V)^n-l,  where  c^res(A ,A  ' )  .  Since  A(x)=anT^=1(x-^) , 


we  have 


,(x)*aX-i 


i/j 


Evaluating  (30)  at  ,  we  obtain 

A  '  (^ )  =anT1i=2  "<^i ) ' 

Hence,  evaluating  (29)  at  x^  and  using  (31): 

{V^^l^lv^c. 


(30) 


(31) 


(32) 


By  [10],  Section  28,  c=anD,  where  D  is  the  discriminant  of  A,  a  non¬ 
zero  Gaussian  integer.  Hence  V(?6^)^0  and  by  (32), 


sep  (Aj-D/V^)^  U  . 

i—  I  i 


| A  '  |  <n  | A [  so 


Ivl^d2"-2 


(33) 


(3*0 


by  Theorem  2.  Since  deg(v)sn-l  and  w  <d , 


o 


W-  .  (35) 

From  (33)  and  (35),  using  |D|al  and  | -q^ | <2d , 

sep(A)22  n  d  (36) 

The  proof  is  completed  by  observing  that  ns2 .  I 

In  order  to  obtain  the  third  root  separation  bound,  we  construct 
a  Gaussian  polynomial  B*  whose  roots  are  all  the  differences 
with  i^ j .  The  idea  of  constructing  B*  as  a  resultant  was  suggested 
by  some  current  research  of  R.  Loos,  [  5] •  After  obtaining  upper 
bounds  for  the  coefficients  of  B*,  we  will  apply  the  following  theorem 
to  obtain  a  lower  bound  for  the  roots  of  B*,  and  hence  for  sep(A) . 

Theorem  6.  Let  A(x)=^_  ^a^x^  be  a  complex  polynomial  of  degree 
n>0,  with  0.  If  &  is  any  root  of  A,  then 

lal^in!^  U^aJ1^  (37) 

a .  a^Q 

1 

Proof.  Let  A*(x)=xnAfx  ^)=^?  „a  ,x^.  A*  is  a  polynomial  of 
-  *_i=0  n-i  r 

degreenwhose  roots  are  the  reciprocals  of  the  roots  of  A,  for 
A*(x)=anxV=1(x_1-  ai)-anrJ.l(l-0rix)-£anTTj-1(-£vi)]^.1)(x-  a.  )}  = 

an(ao/an^TTi=l^x-  ffi  1)=a0T7i=l^X‘°'i  ^  Hence  A^(o'"1)=0  and  from  [!|  J, 
Exercise  4.6.2.20,  we  have 


\a~1\<2  maxUi^|ai/a()|1 \  (30 

from  which  (37)  is  immediate.  M  (39) 

Theorem  7<  ('’ollins,  1973)  Let  A  be  a  Gaussian  poly¬ 
nomial  of  degree  n^2  with  only  simple  roots  and  d= | A | .  Then 

sep(A)^(esr. ’/^d)  n,  (40) 


where  e  is  the  base  of  the  natural  logarithm. 


1 


Proof.  Let  B(x)  be  the  resultant  of  A  y)  and  A  x+y  .  If  the 
coefficients  and  roots  of  A  are  given  by  (1),  then, 

A(x+y)=anT^=1(y-(0j-x)) .  (U) 

Expressing  the  resultant  B(x)  as  a  symmetric  function  of  the  roots  of 
A(y)  and  A(x+y)  by  the  theorem  of  van  der  Waerden  10] ,  Section  2  ), 


2n 


B(x)=an 

Since  o,i=0'j  if  and  onlv  if  i=j ,  B(x)=xnB(x),  where 

B(x)=a^nrri^j(x-(0i-aj)),  (f.j) 

is  a  polynomial  of  degree  n(n-l)  with  B(0)aQ.  Also,  ('*5)  can  be  written 
in  the  form 


B(x)=a^nTTi<j(x2-(ai-aj)'J) ,  0^) 

so  that  if  B(x)=J?_n  ^b.x*  then  1^=0  for  i  odd. 

Expanding  A(x+y)  in  a  Taylor  series, 

A(x+y)=2=0{A(;iVy)/i:}x1,  0*5) 

where  A^  is  the  ir>1  derivative  of  A.  Let 

A*(x,y)  =  rA(x4y)-Afy)l/x=^?=1{A(ii(y)/il)x1  1 .  (H6) 

Let  M  be  the  Sylvester  matrix  of  A(y)  and  A(x+y).  If  we  subtract  the 

i^  row  of  M  from  the  (n+i)  row  and  then  divide  the  latter  by  x,  for 

l^i^n  we  obtain  a  matrix  M  such  that  det(M)=xndet(M) .  The  first  column 

of  M  contains  a  in  the  first  row  and  zeros  elsewhere.  Hence  det(M)= 
n  ' 

a^det(M*),  where  M*  results  from  M  upon  deletion  of  its  first  row  and 
column.  But  M*  is  the  Sylvester  matrix  <>f  A(y)  a  :d  A>  x,y),  so 

Bfx)=anB*(x)  (!•?) 

where  B*(x)  is  the  resultant  of  A  y)  and  A'^Xjy). 
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We  mw  proceed  to  obtain  bounds  for  the  coefficients  of  B‘ .  Let 

\*k+l ,  t  .  .  ,  i-1 


A^k(x,y)-^it1{A  i(y)/iHx 


('*0 


so  that  A *  is  the  result  of  deleting  from  A¥  all  terms  of  degree 
k 

k+1  or  greater  in  x.  Since  k¥  and  A£  are  both  of  degree  n-1  in  y, 

«  I  i 

Bv(x)=B> (x)  (modulo  x  )  for  k20.  Hence  the  coefficients  of  x  in 
B*(x)  and  B£(x)  arc  identical,  and  if  B>*  x  )=^_  b^x  then 


Now  | A  ^(y)/il|^  s0»  by  0*£)> 

lA*(x,y)  U2^\(i)dienk+ld* 

By  Theorem  2  and  (50), 

i du i  n  (k+l)n  2n-l 
|B*|€e  n  d 

By  (U9 ;  and  (51),  together  with  'b*l>l, 

for  k>l.  Since  b*=0  for  i  odd,  by  Theorem  6, 

i(el/?n5/2d)‘n, 


0*9) 


(50) 


(51) 


(52) 


(55) 


completing  the  proof.  I 

The  computing  time  of  an  algorithm,  e.g.  )] ,  for  isolating  the 
zeros  of  a  Gaussian  polynomial  A  is  dominated  in  the  sense  of  by 

a  polynomial  function  of  n=deg  A),  log  d  where  d=|A|,  and  -log  sep  A  . 
If  denotes  codominance  of  functions  as  in  [2]  and  if  n,d  , 
H(n,d)  and  C?(n,d)  are  the  bounds  on  sep  A)  given  by  Theorems  !* ,  and 


7,  then  we  have 


2 

-log  C^(n,d)~ji  log  d. 


(5!*) 


LJ 
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whereas 

-log  H''n,d)~-log  C  (ii,d),jo  log  nd .  (55) 

In  this  *>ense  the  last  two  bounds  are  equivalent. 

2 

When  n=2,  sep^A)  can  be  given  explicitly.  If  A  x)=ax  +bx+c  has  two 
distinct  roots,  then 

sep(A)  =  |b  -  ac|^^/|a|.  (56) 

Also,  by  Theorem  H,  s ’p(A )>l/2d .  l.et  a=k,  b=2k-l  and  c=k-l  with  k>l . 
Then  d=  | A^=U k-2  and  sep(A)*l/k<3;/(l*k-2)*4/d . 

Define 

L(n,d'>=minfsep(A)  :deg(A)=n&l A  I  <d]  .  (57) 

Then,  we  have  just  shown, 

L(L',d')~xfl.  (58) 

It  does  not  seem  unreasonable  to  ask  for  an  explicit  relation  such 
as  (58)  for  L(3,d),  but  we  have  thus  far  not  succeeded  with  this 
apparently  simple  problem.  We  know  only,  by  Theorem  57  and  some 
obvious  examples,  that 

d_5!SL(3,d)^d'1,  (59) 

where  "<"  is  the  dominance  relation. 
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